Theory of anisotropic exchange in laterally coupled quantum dots 



Fabio Baruffa^, Peter Stano^'^ and Jaroslav Fabian-'^ 
^Institute for Theoretical Physics, University of Regensburg, 93040 Regensburg, Germany 

^Institute of Physics, Slovak Academy of Sciences, 84-511 Bratislava, Slovak Republic 
^Physics Department, University of Arizona, 1118 E 4**" Street, Tucson, AZ 85721, USA 

The effects of spin-orbit coupling on the two-electron spectra in lateral coupled quantum dots are 
investigated analytically and numerically. It is demonstrated that in the absence of magnetic field 
the exchange interaction is practically unaffected by spin-orbit coupling, for any interdot coupling, 
boosting prospects for spin-based quantum computing. The anisotropic exchange appears at finite 
magnetic fields. A numerically accurate effective spin Hamiltonian for modeling spin-orbit-induced 
two-electron spin dynamics in the presence of magnetic field is proposed. 

PACS numbers: 71.70.Gm, 71.70.Ej, 73. 21. La, 75.30.Et 



The electron spins in quantum dots are natural and 
viable qubits for quantum computing, [H as evidenced by 
the impressive recent experimental progress 0, Q in spin 
detection and spin relaxation, as well as in coherent 
spin manipulation. 0] In coupled dots, the two-qubit 
quantum gates are realized by manipulating the exchange 
coupling which originates in the Coulomb interaction and 
the Pauli principle. [ll,@| How is the exchange modified by 
the presence of the spin-orbit coupling? In general, the 
usual (isotropic) exchange changes its magnitude while 
a new, functionally different form of exchange, called 
anisotropic, appears, breaking the spin-rotational sym- 
metry. Such changes are a nuisance from the perspective 
of the error correction, [91 although the anisotropic ex- 
change could also induce quantum gating, l^, Hi 

The anisotropic exchange of coupled localized electrons 
has a convoluted history 12 - 13| • The question boils down 
to determining the leading order in which the spin-orbit 
coupling affects both the isotropic and anisotropic ex- 
change. At zero magnetic field, the second order was 
suggested, [3| with later revisions showing the effects are 
absent in the second order, [l^, [3] The analytical com- 
plexities make a numerical analysis particularly useful. 

Here we perform numerically exact calculations of the 
isotropic and anisotropic exchange in realistic GaAs cou- 
pled quantum dots in the presence of both the Dres- 
selhaus and Bychkov-Rashba spin-orbit interactions. 21 
The numerics allows us to make authoritative statements 
about the exchange. We establish that in zero magnetic 
field the second-order spin-orbit effects are absent at all 
interdot couplings. Neither is the isotropic exchange af- 
fected, nor is the anisotropic exchange present. At finite 
magnetic fields the anisotropic coupling appears. We de- 
rive a spin-exchange Hamiltonian describing this behav- 
ior, generalizing the existing descriptions; we do not rely 
on weak coupling approximations such as the Heitler- 
London one. The model is proven highly accurate by 
comparison with our numerics and we propose it as a re- 
alistic effective model for the two-spin dynamics in cou- 
pled quantum dots. 

Our microscopic description is the single band effec- 



tive mass envelope function approximation; we neglect 
multiband effects. 22, [2^ We consider a two electron 
double dot whose lateral confinement is defined electro- 
statically by metallic gates on the top of a semiconduc- 
tor heterostructure. The heterostructure, grown along 
[001] direction, provides strong perpendicular confine- 
ment, such that electrons are strictly two dimensional, 
with the Hamiltonian (subscript i labels the electrons) 



H=Y.{T, + V, + Hz,^ + -ffscO + He 



(1) 



The single electron terms are the kinetic energy, model 
confinement potential, and the Zeeman term, 

T = PV2to= (-i7iV + eA)V2m, (2) 
V = {l/2)mu'^[u):m{{x-df,{x + df} + y\ (3) 
Hz = {g/2){eh/2me)B ■ cr = fiB ■ CT, (4) 

and spin-orbit interactions — linear and cubic Dressel- 
haus, and Bychkov-Rashbafil) . 
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which we lump together as Hgo — mv ■ cr. The posi- 
tion r and momentum P vectors are two dimensional 
(in-plane); m/rrie is the effective/electron mass, e is the 
proton charge, A — Bz{~y, x)/2 is the in-plane vector 
potential to magnetic field B = {Bx, By, Bz), g is the 
electron g-factor, cr are Pauli matrices, and fx is the renor- 
malized magnetic moment. The double dot confinement 
is modeled by two equal single dots displaced along [100] 
by ±d, each with a harmonic potential with confinement 
energy huj. The spin-orbit interactions are parametrized 
by the bulk material constant 7^ and the heterostruc- 
ture dependent spin-orbit lengths 4r, Id- Finally, the 
Coulomb interaction is He = (e^/47re)|ri — r2|^^, with 
the dielectric constant e. 

The numerical results are obtained by exact diagonal- 
ization (configuration interaction method). The two elec- 
tron Hamiltonian is diagonalized in the basis of Slater 
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FIG. 1: Calculated double dot spectrum as a function of the 
interdot distance/tunneling energy. Spin is not considered 
and the magnetic field is zero. Solid lines show the two elec- 
tron energies. The two lowest states are explicitly labeled, 
split by the isotropic exchange J and displaced from the near- 
est higher excited state by A. For comparison, the two lowest 
single electron states are shown (dashed), split by twice the 
tunneling energy T. State spatial symmetry is denoted by 
darker (symmetric) and lighter (antisymmetric) lines. 



determinants constructed from numerical single electron 
states in the double dot potential. Typically we use 21 
single electron states, resulting in the relative error for 
energies of order 10~^. We use material parameters of 
GaAs: m = 0.067me, g = -0.44, jc = 27.5 meVA^, a 
typical single dot confinement energy huj = 1.1 meV, and 
spin-orbit lengths Id = 1.26 /Ltm and = 1.72 fim from a 



fit to a spin relaxation experiment. 24l. |25| 

Let us first neglect the spin and look at the spectrum 
in zero magnetic field as a function of the interdot dis- 
tance (2(i)/tunneling energy, Fig. [1] At d = our model 
describes a single dot. The interdot coupling gets weaker 
as one moves to the right; both the isotropic exchange 
J and the tunneling energy T decay exponentially. The 
symmetry of the confinement potential assures the elec- 
tron wavefunctions are symmetric or antisymmetric upon 
inversion. The two lowest states, 5'±, are separated from 
the higher excited states by an appreciable gap A, what 
justifies the restriction to the two lowest orbital wave- 
functions for the spin qubit pair at a weak coupling. Our 
further derivations are based on the observation 



/l/2*± = i^-i. 



(8) 



where If(x,y) = /(— — y) is the inversion operator 
and Pfig-i = /2<?i is the particle exchange operator. 
Functions in the Heitler-London approximation fulfill 
Eq. IH). However, unlike Heitler-London, Eq. ([8]) is valid 
generally in symmetric double dots, as we learn from nu- 
merics (we saw it valid in all cases we studied). 



Let us reinstate the spin. The restricted two qubit 
subspace amounts to the following four states {S stands 
for singlet, T for triplet), 

.,4 - {^+S, *_r+, ^-^To, (9) 

Within this basis, the system is described by a 4 by 4 
Hamiltonian with matrix elements (i?4)ij = {^i\H\^j). 
Without spin-orbit interactions, this Hamiltonian is di- 
agonal, with the singlet and triplets split by the isotropic 
exchange J,[ll,[l| and the triplets split by the Zeeman en- 
ergy ^B. It is customary to refer only to the spinor part 
of the basis states, using the sigma matrices, resulting in 
the isotropic exchange Hamiltonian, 



i/iso = (J/4)cri • 0-2 + • (cTi + (T2) 



(10) 



A naive approach to include the spin-orbit interaction 
is to consider it within the basis of Eq. This gives 
the Hamiltonian H' = H,. 



Kniso-^ where 



H' 



((Ti - (T2) + b' • (cTi X (T2) 



(11) 



with the six real parameters given by spin-orbit vectors 

a' = Rc(«'+|wi|1'_), b' = Im(«'+|wi|1'_). (12) 

The form of the Hamiltonian follows solely from the in- 
version symmetry /w = — w and Eq. ([8]). The spin-orbit 
coupling appears in the first order. 

The Hamiltonian H[,^ fares badly with numerics. Fig- 
ure [2] shows the energy shifts caused by the spin-orbit 
coupling for selected states, at different interdot cou- 
plings and perpendicular magnetic fields. The model is 
completely off even though we use numerical wavefunc- 
tions ± in Eq. (jl2p without further approximations. 

To improve the analytical model, we remove the 
linear spin-orbit terms from the Hamiltonian using 
transformationjiO, S 

U = exp [-(i/2)ni • (Ti - (i/2)n2 • (T2] , (13) 

where n = {x/U - y/hr, x/kr - y/k, 0). 

Up to the second order in small quantities (the spin- 
orbit and Zeeman interactions), the transformed Hamil- 
tonian H = UHU^ is the same as the original, Eq. ([T}, 
except for the linear spin-orbit interactions: 

Hso = -(/iB X n) • cr + {K^/n)L,a, - K+, (14) 

where K± = {fi^ /Amll) ± (h^/Amll^). In the unitarily 
transformed basis, we again restrict the Hilbert space to 
the lowest four states, getting the effective Hamiltonian 



-ffcx = (^/4)cri • 0-2 + fi{B + B,o) ■ (cTi + 0-2) 
+ a • (cTi - 0-2) + b • (cTi X 0-2) - 2K+. 



(15) 



The operational form is the same as for H^^. The qual- 
itative difference is in the way the spin-orbit enters the 
parameters. First, a contribution to the Zeeman term. 



(16) 
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FIG. 2: The spin-orbit induced energy shift as a function of 
the interdot distance (left) and perpendicular magnetic field 
(right), a) Singlet in zero magnetic field, c) singlet at 1 Tesla 
field, b) and d) singlet and triplet T+ at the interdot distance 
55 nm corresponding to the zero field isotropic exchange of 
1 /ieV. The exchange models H[,^ (dashed) and _ffex (dot- 
dashed) are compared with the numerics (solid). 



appears due to the inversion symmetric part of Eq. ([H]) . 
Second, the spin-orbit vectors are hnearly proportional 
to both the spin-orbit coupling and magnetic field, 

a = -/iB X Rc(1'+|ni|^'_), (17a) 
b = -/iB X Im(*+|ni|*_). (17b) 

The effective model and the exact data agree very well 
for all interdot couplings, as seen in Fig. O 

At zero magnetic field, only the first and the last term 
in Eq. (fTS]) survive. This is the result of Ref [2^, where 
primed operators were used to refer to the fact that the 
Hamiltonian H^-x refers to the transformed basis, {U^i}. 
Note that if a basis separable in orbital and spin part is 
required, undoing U necessarily yields the original Hamil- 
tonian Eq. (U]), and the restriction to the four lowest 
states gives H'^^. Replacing the coordinates {x,y) by 
mean values (id, 0)[l3| visualizes the Hamiltonian H^x 
as an interaction through rotated sigma matrices, but 
this is just an approximation, valid if d, I so ^ ^o- 

One of our main numerical results is establishing the 
validity of the Hamiltonian in Eq. (ITSI) for i? = 0, con- 
firming recent analytic predictions and extending their 
applicability beyond the weak coupling limit. In the 
transformed basis, the spin-orbit interactions do not lead 
to any anisotropic exchange, nor do they modify the 
isotropic one. In fact, this result could have been antic- 
ipated from its single-electron analog: at zero magnetic 
field there is no spin-orbit contribution to the tunneling 
energy, [28| going opposite to the intuitive notion of the 
spin-orbit coupling induced coherent spin rotation and 
spin-fiip tunneling amplitudes. Figure summarizes 
this case, with the isotropic exchange as the only nonzero 
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FIG. 3: a) The isotropic and anisotropic exchange as func- 
tions of the interdot distance at zero magnetic field, b) The 
isotropic exchange J, anisotropic exchange c/c', the Zeeman 
splitting /i_B, and its spin-orbit part fiBso at perpendicular 
magnetic field of 1 Tesla. c) Schematics of the exchange-split 
four lowest states for the three models, Hiao, H^^, and //cx, 
which include the spin-orbit coupling in no, first, and second 
order, respectively, at zero magnetic field (top). The latter 
two models are compared in perpendicular and in-plane mag- 
netic fields as well. The eigenenergies are indicated by the 
solid lines. The dashed lines show which states are coupled 
by the spin-orbit coupling. The arrows indicate the redistri- 
bution of the couplings as the in-plane field direction changes 
with respect to the crystallographic axes (see the main text). 



parameter of model Hex- In contrast, model H^^ predicts 
a finite anisotropic exchange. [3(1 

From the concept of dressed qubits[29i| it follows that 
the main consequence of the spin-orbit interaction, the 
transformation U of the basis, is not a nuisance for quan- 
tum computation. We expect this property to hold also 
for a qubit array, since the electrons are at fixed positions 
without the possibility of a long distance tunneling. How- 
ever, a rigorous analysis of this point is beyond the scope 
of this article. If electrons are allowed to move, U results 
in the spin relaxation, [soj 

Figure [31d shows model parameters in 1 Tesla perpen- 
dicular magnetic field. The isotropic exchange again 
decays exponentially. As it becomes smaller than the 
Zeeman energy, the singlet state anticrosses one of the 
polarized triplets (seen as cusps on Fig. [5]). Here it is 
T+, due to the negative sign of both the isotropic ex- 
change and the g-factor. Because the Zeeman energy 
always dominates the spin-dependent terms and the sin- 
glet and triplet Tq are never coupled (see below), the 
anisotropic exchange influences the energy in the sec- 
ond order. [l2| Note the difference in the strengths. In 
i/px the anisotropic exchange falls off exponentially, while 
Hex predicts non-exponential behavior, resulting in spin- 
orbit effects larger by orders of magnitude. The effective 
magnetic field Bso is always much smaller than the real 
magnetic field and can be neglected in most cases. 
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Figure[3]: compares analytical models. In zero field and 
no spin-orbit interactions, the isotropic exchange Hamil- 
tonian iJjso describes the system. Including the spin- 
orbit coupling in the first order, H^^, gives a nonzero 
coupling between the singlet and triplet Tq. Going to the 
second order, the effective model Hc^ shows there are no 
spin-orbit effects (other than the basis redefinition). 

The Zeeman interaction splits the three triplets in a fi- 
nite magnetic field. Both H^^ and iJex predict the same 
type of coupling in a perpendicular field, between the 
singlet and the two polarized triplets. Interestingly, in 
in-plane fields the two models differ qualitatively. In 
ifgx the spin-orbit vectors are fixed in the plane. Rota- 
tion of the magnetic field "redistributes" the couplings 
among the triplets. (This anisotropy with respect to 
the crystallographic axis is due to the C2v symmetry of 
the two-dimensional electron gas in GaAs, imprinted in 
the Bychkov-Rashba and Dresselhaus interactions. [2]|) 
In contrast, the spin-orbit vectors of Hex are always per- 
pendicular to the magnetic field. Remarkably, aligning 
the magnetic field along a special direction (here we al- 
low an arbitrary positioned dot, with S the angle between 
the main dot axis and the crystallographic x axis). 



[Ibr — Id tan S, Id — hr tan 6, 0] , 



(18) 



all the spin- orbit effects disappear once again, as if B 
were zero. (An analogous angle was reported for a sin- 
gle dot in Ref. [3l|). This has strong implications for 
the spin-orbit induced singlet-triplet relaxation. Indeed, 
S' ■f-!> To transitions are ineffective at any magnetic field, 
as these two states are never coupled in our model. Sec- 
ond, 5 o T± transitions will show strong (orders of mag- 
nitude) anisotropy with respect to the field direction, 
reaching minimum at the direction given by Eq. (|18p. 
This prediction is straightforwardly testable in experi- 
ments on two electron spin relaxation. 

Our derivation was based on the inversion symmetry 
of the potential only. What are the limits of our model? 
We neglected third order terms in Hgo and, restricting 
the Hilbert space, corrections from higher excited orbital 
states. (Among the latter is the non-exponential spin- 
spin coupling |12|). Compared to the second order terms 
we keep, these are smaller by (at least) d/lso and c/A, 
respectively. [35] Apart from the analytical estimates, the 
numerics, which includes all terms, assures us that both 
of these are negligible. Based on numerics we also con- 
clude our analytical model stays quantitatively faithful 
even at the strong coupling limit, where A — 0. More 
involved is the infiuence of the cubic Dresselhaus term, 
which is not removed by the unitary transformation. 
This term is the main source for the discrepancy of the 
model and the numerical data in finite fields. Most im- 
portantly, it does not change our results for B = 0. 

Concluding, we studied the effects of spin-orbit cou- 
pling on the exchange in lateral coupled GaAs quantum 



dots. We derive and support by precise numerics an effec- 
tive Hamiltonian for two spin qubits, generalizing the ex- 
isting models. The effective anisotropic exchange model 
should be useful in precise analysis of the physical realiza- 
tions of quantum computing schemes based on quantum 
dot spin qubits, as well as in the physics of electron spins 
in quantum dots in general. Our analysis should also 
improve the current understanding of the singlet-triplet 
spin relaxation 32h35|. 
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